Mott Transition, Antiferromagnetism, and d-wave Superconductivity in 
Two-Dimensional Organic Conductors 
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We study the Mott transition, antiferromagnetism and superconductivity in layered organic 
conductors using Cellular Dynamical Mean Field Theory for the frustrated Hubbard model. A 
d-wave superconducting phase appears between an antiferromagnetic insulator and a metal for 
t'/t = 0.3 — 0.7, or between a nonmagnetic Mott insulator (spin liquid) and a metal for t' /t > 0.8, 
in agreement with experiments on layered organic conductors including k-(ET)2Cu2(CN)3. These 
phases are separated by a strong first order transition. The phase diagram gives much insight into 
the mechanism for d-wave superconductivity. Two predictions are made. 

PACS numbers: 71.10.Fd, 71.27.-fa, 71.30.-fh, 71.10.-w 
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Strong electronic correlations lead to fascinating phe- 
nomena such as high-temperature superconductivity and 
metal-insulator (Mott) transitions. In that context, lay- 
ered organic conductors k-(BEDT-TTF)2X (X denotes 
an anion) play a special role. They share many features 
with high-temperature superconductors (HTSC) 0, 13 j 
such as the existence of d-wave superconductivity and 
antiferromagnetism. But they have an even richer phase 
diagram since they also display a Mott transition and 
possibly a T = spin liquid phase. Microscopically, they 
are described by the two-dimensional one-band Hubbard 
model, just as HTSC, albeit on an anisotropic triangular 
lattice instead of a square lattice. (2| The various phases 
of HTSC are explored essentially by changing doping. 
The phases of organic conductors on the other hand are 
controlled by pressure and by frustration, two variables 
over which one has very little control in the HTSC. In this 
paper, we will see that we can capitalize on recent the- 
oretical progress on the theory of HTSC to understand 
organic conductors. In return, the agreement that we 
find between theory and experiment for these compounds 
leads to much insight into the origin of d-wave supercon- 
ductivity in the one-band Hubbard model in general. 

Experimentally, metallic, paramagnetic insulating 
(spin liquid SL), antiferromagnetic (AF), and unconven- 
tional superconducting (SC) phases are all found, for ex- 
ample, in the pressure vs temperature phase diagram 
of X=Cu[N(CN)2]Cl Q. The broken-symmetry states 
overlap through a first order boundary that merges with 
the first order line of the Mott metal-insulator transi- 
tion in the normal state. Various experiments suggest 
that the superconductivity has line nodes 5] that are 
strongly suggestive of d-wave character. Changing the 
anion modifies the frustration in the lattice. The anion 
(X=Cu2(CN)3) Q, which corresponds to large frustra- 
tion, has attracted considerable attention because of the 
transition from d-wave SC to a possible SL state where 
no magnetic long-range order was found down to a very 
low temperature (23mK), in contrast to other layered or- 
ganic materials. 

The intriguing behaviors observed in layered organic 



conductors have prompted theoretical studies using var- 
ious analytical and numerical methods. The existence 
of magnetic, metallic and spin liquid phases, but not d- 
wave superconductivity, has been explored using a path- 
integral renormalization group method for the frustrated 
Hubbard model by Morita et al. • Signatures of a finite 
temperature Mott critical point in agreement with exper- 
imental studies of K-organics and with single site DMFT 
have been found with Cellular Dynamical Mean Field 
Theory (CDMFT) S\ in conjunction with a quantum 
Monte Carlo (QMCJ method, by ParcoUet et al. ^]. At 
T = 0, using a variational Monte Carlo (VMC) method, 
Liu et al. 'lO'| showed an unconventional SC ground state, 
sandwiched between a conventional metal at weak cou- 
pling and a spin liquid at large coupling, but antiferro- 
magnetism was not considered. Similarly, with a U(l) 
gauge theory in the slave-rotor representation, Lee et 
al. Illl found a first order transition from a supercon- 
ductor to a spin liquid and Powell et al. 12] showed a 
first order transition from a Mott insulator to a d-wave 
superconductor using a Gutzwiller projection method for 
the Hubbard-Heisenberg model. Can et al. [l3j found a 
Gossamer superconductor at small U and an AF insula- 
tor at large U, separated by a first order transition. Only 
the variational study of Watanabe et al. jl4] considered 
the possibility of all phases, metallic, d-wave SC, AF and 
SL. However, they did not find AF where it is observed 
Q experimentally. 

The layered organic conductors are described by the 
frustrated two-dimensional Hubbard model Il5l 
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where cj^ (cia) are creation (annihilation) operators for 
electrons of spin a, rii^ = cl^Cia is the density of 
a spin electrons. The hopping amplitude tij deter- 
mines the anisotropic bare dispersion e^ ~ — 2i(cos k^ + 
cos ky) — 2t' cos{kx + ky). (The HTSC have an additional 
2t' cos{kx — ky) in their dispersion). U is the on-site re- 
pulsive interaction and /i is the chemical potential con- 
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FIG. 1: Phase diagram of the frustrated two-dimensional 
Hubbard model at zero temperature in the C//f vs t'/f plane. 
AF, SL, SC and M denote antiferromagnetic, spin liquid 
(nonmagnetic insulator), d-wave superconducting and metal- 
lic phases. 



trolling the electron density. Physically relevant param- 
eters for layered organic conductors are U jt = 5 — 10 at 
half-filling (n = 1) and i'/i = 0.5 - 1.1. 

Our results are summarized by the \J jt vs t' jt phase 
diagram in Fig. ^ A d-wave SC phase appears in a 
relatively narrow region between metallic (M) and AF 
phases, or between metallic (M) and spin liquid phases 
(SL). The phase transition lines are all first order, ex- 
cept for the SC-metal transition that is also first order 
but eventually turns second order (dashed line) for suf- 
ficiently large t' jt (see caveat flq). Three phases, M- 
AF-SC, meet at the triple point near i'/i = 0.25 and 
V jt — 2.75, satisfying the Gibbs phase rule. There is a 
jump for all order parameters so that this point does not 
exhibit SO (5) symmetry [l3|. The SL phase (paramag- 
netic Mott insulator) appears only for large t' jt (large 
frustration) and large U jt (strong correlations). If we do 
not allow for broken symmetry, we obtain a normal state 
that has a first order Mott metal-insulator transition that 
coincides with the first order line that separates SC from 
either AF or SL. That first order line however ends at a 
critical point for i'/t close to 0.5. 

Predictions : Based on the calculated phase diagram, 
we predict a class of new materials (with t' /t ^ 0.8 — 0.9) 
which would undergo a series of phase transitions from 
an AF insulator to a paramagnetic Mott insulator (spin 
liquid) to a d-wave superconductor to a metal with in- 
creasing pressure. Also, we predict that if AF and SC 
are destroyed with an external symmetry-restoring per- 
turbation (magnetic field) , the first-order transition that 
separates them should remain as a first-order Mott metal- 
insulator transition. In the following, we describe the 
method and results that lead to the above phase dia- 
gram. We then discuss the remarkable agreement with 



existing experiments. 

Method : We use Cellular Dynamical Mean-Field The- 
ory (CDMFT) 8], a cluster approach that allows one 
to extend Dynamical Mean Field Theory (DMFT) to in- 
clude momentum dependence of the self-energy. CDMFT 
has been benchmarked and is accurate even in one dimen- 
sion [13, [l3l ■ The infinite lattice is tiled with identical 
clusters of size TVc, and the degrees of freedom in the 
cluster are treated exactly while the remaining ones are 
replaced by a bath of non-interacting electrons that is 
determined self-consistently. To solve the quantum clus- 
ter embedded in an effective SC medium, we consider a 
cluster-bath Hamiltonian of the form 123, EE 
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Here the indices ^,v — 1 
cluster, and c„a and a\ 



, Nc label sites within the 
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ter and the bath, respectively, t 



annihilate electrons on the clus- 
^1, is the hopping matrix 



within the cluster, and e'^^ are the bath energies and 
^m^CT ^I'S the bath-cluster hybridization matrices. Am„ 
represents the amplitude of SC correlations on the bath 
with a given gap symmetry. Because superconductivity 
and antiferromagnetism are allowed to compete on equal 
footing, e^g. and F^^g. carry a spin variable a explic- 
itly. In the present study we used A^c = 4 sites for the 
cluster and A^ = 8 sites for the bath with m = 1, ...,4, 
a = 1,2. In the normal state, because of symmetry, 
there are 8 independent parameters, while there are 18 
when AF and SC are allowed to compete. To deal with 
superconductivity, the Nambu spinor representation is 
used for the cluster operators so that the Weiss field, 
the cluster Green's function and self-energy constructed 
from these operators are 8x8 matrices. The exact diag- 
onalization method 22] is used to solve the cluster-bath 
Hamiltonian Eq. [3 at zero temperature, which has the 
advantage of computing dynamical quantities directly in 
real frequency and of treating the large U regime without 
difficulty. Although the present study focuses on a 2 x 2 
cluster with additional 8 bath sites, we expect our re- 
sults to be robust with respect to an increase in the clus- 
ter size. This was verified by our recent low (but finite) 
temperature CDMFT-I-QMC calculations [23| where at 
intermediate to strong coupling a 2 x 2 cluster accounts 
for more than 95% of the correlation effect of the infi- 
nite size cluster in the single-particle spectrum. Recent 
Variational Cluster Perturbation Theory (VCPT) calcu- 
lations 24] for the same Hamiltonian also confirmed that 
results on a 2 X 2 cluster are quantitatively similar to 
those on larger clusters. 

Normal state Mott transition: We first present the evi- 
dence of a first order line of metal-to-insulator transition 
in the normal state. That line also turns out to coincide 
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FIG. 2: (a) double occupancy D = (n|n|) as a function of 
U/t for t'/t = 0.7 (solid curve) and t' /t = 1.0 (dashed curve), 
(b) critical Uc/t as a function of t'/t in the normal state. 
The diamond symbols in (a) are obtained with the insulating 
solution as an initial guess. 



with the phase boundary between AF and SC phases for 
t'/t > 0.5. The present results obtained in the normal 
state at r = would be relevant when the temperature 
is slightly higher than Tc and T/v, or in regions where 
broken-symmetries are destroyed at zero temperature by 
an external symmetry restoring perturbation. The first 
plot in Fig.|2a) shows the double occupancy (nini) (cir- 
cles) for two different t'/t. For t'/t — 0.7 it jumps dis- 
continuously to a lower value near U/t = 7.0, while for 
t'/t = 1.0 the first order transition occurs near U/t = 8.5 
which is comparable to the value obtained by ParcoUet 
et al. Ig| in CDMFT-hQMC and Watanabe et al. |jj| in 
VMC. To show hysteresis associated with a first order 
transition, we also calculated the double occupancy (di- 
amonds in Fig. EJa)) with the insulating solution as an 
initial guess. In this case the first order transition occurs 
for a value of U that is typically smaller by 0.5 — 0.75i. 
A very similar hysteresis is found for all the first order 
transitions shown in this paper, but for the rest of the 
paper we show results obtained only with the metallic 
solution as an initial guess. The critical Uc shown in 
Fig. Efb) decreases monotonically with decreasing t'/t 
until the first order transition disappears near t'/t — 0.5, 
close to t'/t < 0.4 found by Watanabe et al. Q. 

While various studies [lH El El El El suggested 
the existence of a first order transition for the frustrated 
Hubbard model in static quantities, here we present 
the strongest evidence of a first order transition from a 
metal to an insulator by presenting dynamical quantities 
A{k, oj). For U/t = 6.9 slightly smaller than critical Uc/t 
for t'/t — 0.7, the spectral function has a sharp quasipar- 
ticle peak near the Fermi wave vectors (Fig.|3Ia)) consis- 
tent with a Fermi liquid picture. It is more evident from 
the fact that A{k,uj) becomes sharper as k approaches 
the Fermi surface. When U/t is increased only by a tiny 
fraction (less than 1%), the first order metal-to-insulator 
transition manifests itself as the massive reshuffling of 
the spectral weight. An insulating gap is present in the 



FIG. 3: Spectral function A{k,uj) along some symmetry di- 
rections for (a) U/t = 6.9 and (b) U/t = 6.95 for t'/t = 0.7 in 
the normal state. 
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FIG. 4: (Color online) Order parameters of the frustrated 
two-dimensional Hubbard model as a function of pressure 
(modeled as t/U) at zero temperature for (a) t'/t — 0.7 and 
(b) t'/t — 1.0. d (circles) and m (diamonds) denote the d-wave 
SC and AF order parameters, respectively, m is multiplied 
by 0.1 to fit in the plots. 



whole Brillouin zone, reminiscent of the spectral function 
in the two-dimensional Hubbard model at half-filling |20| 
where the low energy bands inside the Hubbard bands 
are caused by short-range spin correlations. 

Broken symmetry states: Next we study the existence 
of broken symmetry states in the frustrated model and 
the role of the first order transition in those states. 
Because of possible coexistence of AF and SC phases, 
the two order parameters are treated on equal footing. 
For a SC phase two symmetries are considered dx2-y2 
and d^y. For all the parameters studied d^2_„2 sym- 
metry always dominates over d^y (see caveat '2^). For 
t'/t = 0.7 (Fig. Ufa)) AF and SC phases abruptly termi- 
nate near U/t — 7, the same Uc as in the normal state 
where a metal-to-insulator transition occurs. We have 
not found evidence of homogeneous states with non-zero 
AF and SC order parameters simultaneously, in contrast 
to the case of the square lattice t — t' — U model for 
HTSC at finite doping, where AF can have some itinerant 
character |2^|23. With enough frustration, t'/t =1.0 
(Fig. EJb)) a direct first order transition from a param- 



agnetic Mott insulator to a superconductor occurs near 
U/t = 8.25 in agreement with other recent theoretical 
works [ill Il2 . [lj|. That paramagnetic Mott insulator 
is commonly refered to as a spin liquid. This result 
is also in agreement with the study [23 of the J — J' 
Heisenberg model in which the Neel order persists up to 
J' /J = 0.6 - 0.7 {t'/t = 0.77 - 0.84). Nevertheless, we 
have to leave open the possibility that magnetically or- 
dered phases different from AF are more stable than the 
SL. 

Contact with experiment: In first approximation, the 
effect of pressure is to decrease the ratio U/t, Q even 
though t'/t should also change slightly [23. We now 
demonstrate that our phase diagram Fig. Q] accounts for 
a surprisingly large number of experimental results if we 
assume that pressure only decreases the U/t ratio at 
fixed t'/t, the various compounds being associated to 
a given t'/t. For t'/t = 0.7 (Fig. g^a)) which is rel- 
evant to K-(ET)2Cu[N(CN)2]Cl, Fig. [2 shows that, as 
observed experimentally Q|: (a) as pressure increases, 
{U/t decreases), one crosses a first order transition be- 
tween AF and SC phase (b) the maximum SC order pa- 
rameter occurs at the phase boundary (c) in the normal 
state there is a first order Mott-insulator transition, at es- 
sentially the same t/U ratio. The d-wave SC in our phase 
diagram however exists in a relatively small region of t/U 
in Fig. 01 In the experimental phase diagram |^J the SC 
region extends far beyond the AF+SC phase boundary 
in the pressure vs temperature plane, but in the absence 
of a precise scale connecting pressure and t/U this can- 
not be taken as a real disagreement. For t'/t — 1.0, close 
to the value for k-(ET)2Cu2(CN)3, the AF phase is not 
stabilized even for large U/t due to too strong frustration 
and the transition is between a SL and SC, (Fig. Elb)) 
consistent with recent experiments by Shimizu et al. ,Q|. 

Another clear overall trend in our results is that 
the maximum value of the SC order parameter de- 
creases monotonically with t'/t (0.0276, 0.0258, 0.0228, 
0.0186, 0.0145 for t'/t = 0.6,0.7,0.8,0.9,1.0). This 
trend is remarkably consistent with experiments where 
compounds with weaker frustration (smaller t'/t) have 
higher T^, for instance, [1 H T^ = 11.6, 10.4, 3.9K for 
for X=Cu[N(CN)2]Br(t7i = 0.68), X=Cu(NCS)2(i7t = 
0.84), X=Cu2(CN)3(i7t = 1.06), respectively. 



Mechanism for d-wave superconductivity: Along the 
AF-SC and SL-SC phase transition lines, the maximum 
value of the SC order parameter increases with U un- 
til U ^ 6t and then it starts to decrease, as does the 
value of J = 4iVC/. In previous studies of HTSC, .2J.,id 
the maximum value of the doping dependent SC order 
parameter was found to scale similarly with interaction 
strength. In the organics, at fixed t'/t the SC order pa- 
rameter always increases with increasing U until it drops 
to zero at the phase boundary. The SC region is broadest 
near t'/t = 0.5 and becomes narrower for both smaller 
and larger t'/t, vanishing near t'/t = 0.25 and t'/t — 1.2, 
respectively. If there is not enough frustration in the 
system, as in the case of near perfect nesting at small 
t'/t, AF long-range order is stabilized, suppressing the 
SC phase. Increasing frustration then helps d-wave SC, 
as suggested for the cuprates |31|. However, too strong 
frustration at large t'/t suppresses even short-range sin- 
glet correlations on which pairing correlations may build 
up so that SC disappears again for t'/t > 1.2. Note that 
the ratio of interaction strength U to hopping integral t 
is of the same order in the organics and in HTSC but 
that the energy scale set by the hopping integrals is an 
order of magnitude smaller in these compounds than in 
the HTSC. 

In conclusion, we obtained the phase diagram for lay- 
ered organic conductors using CDMFT. The calculated 
sequence of phases and the nature of the transitions be- 
tween them are consistent with observations in that class 
of compounds. In addition, the observed maximum value 
of Tc near the Mott transition is consistent with the cal- 
culated maximum SC order parameter. This allows us 
to predict a new class of materials. In conjunction with 
results for HTSC |23, 1^3 1 ^^^ present results also allow 
us to understand more deeply the role of exchange inter- 
actions and of frustration for d-wave superconductivity. 
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